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Higher anisotropic d-wave symmetry in cuprate superconductors 
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We derive a pair potential from tight binding further neighbours attraction that leads to super- 
conducting gap symmetry similar to that of the phenomenological spin fluctuation theory of high 
temperature superconductors (Monthoux, Balatsky, Pines, Phys. Rev. Lett. 67, 3448). We show 
that higher anisotropic d-wave than the simpliest d-wave symmetry is one of the important ingredi- 
ents responsible for higher BCS characteristic ratio. 



PACS Numbers: 74.20Mn, 74.62-c, 74.25DW 

Pairing symmetry of the superconducting energy gap 
in high temperature superconductors still remains an 
open problem after a decade of its discovery. Various ex- 
perimental results which lead to conflicting conclusions 
resulted no concrete concencus to the theory of pairing 
mechanism for high T c superconductors. However, there 
are strong evidences that the pairing state of the cuprate 
superconductors could be d-wave like; experimental ,-ob- 
servations that are sensitive to the phaseEJ and nodesa of 
the gap, reported a sign reversal of the order parameter 
supporting d-wave pairing. On the other hand, a group 
of experiments on the same YBaiCu^Oi (YBCO).mate- 
rial indicate existence of a significant s-componentn; had 
YBCO been a pure d-wave superconductor, it would be 
orthogonal to s-wave state of Pb resulting zero Joseph- 
son supercurrent (while in experiment a well defined c- 
axis current is seen) and there are strong evidences that 
the electron-doped Nd2- x Ce x Cu04 superconductors are 
s-wave typeo. There are also indications, both from theo- 
ries and experiments, that the high T c materials may have 
a mixed pairing symmetry (e.g, d±s or d+is/d xy etc.) ki 
presence of externaLmagnetic field, magnetic impurityS, 
interface effects etcB In addition, there exists important 
clues that indicate pairing state even in the bulk of the 
cuprates and in absence of magnetic field may also have 
a mixed pairing state, with a minor component coex- 
isting with predominant d-wave (i.e., d + e z9 a scenario, 
a = s, d xy )u. The angle resolved photoemission spec- 
troscopy (ARPES) study by Kelley et al, provides strong 
indication that Bi2Sr2CaCu20s+i compound is d-wave 
like in the under- and optimally doped regime whereas 
not a d-wave like in the slightly overdoped but high-T c 
sample. Raman measurements confirmed the unexpected 
behavior of gap symmetry (from predominant d x i_ y i in 
under- or optimally doped to anisotropic s-wave type in 
overdoped) by overdoping Bi2Sr2CaCu20s+5 and almost 
similar phenomena is also found in other high-T c material 
Tl 2 Ba2Cu06+,5El. Now, observation of any s component 
will have stringent contraints on various potential theo- 
ries of high T c superconductors, such as, antiferromag- 
netic spin fluctuation theoryH. Because, antiferromag- 
netic spin fluctuation theory leads to attraction in the 
d-wave channel and pair breaking in the s-wave chan- 
nel (but in a model calculation via spin fluctuation in 
heavy fermion systems by Miyake et alS3 indicates possi- 



bility of anisotropic s-wave pairing as well) . On the other 
hand, pairing mechanisms based on clcctron-phonon in- 
teractions, polarons etc. would be compatible with pure 
d-wave, pure s-wave or an admixture of the twoEU. 

Therefore, it is evident from the above discussion that 
the symmetry of the order parameter and the associated 
mechanism of pairing in high T c cuprates are not at all 
clear, but is essential to have a first step development to- 
wards an understanding to the mystry of pairing mecha- 
nism. The spin fluctuation theory is one of the potential 
theories of high T c superconductivity which can account 
a number of anamolous properties observed in cuprates. 
It is a phenomenological theory with a few small param- 
eters, like the phenomenological form of the spin sus- 
ceptibility x(Q)> f ne magnetic coherence length £, the 
magnon frequency ws p extracted from Nuclear Magnetic 
Resonance (NMR) experiment^. The approximate mo- 
mentum distribution of the superconducting energy gap 
function obtained by the authors of ref.Q is, 



&SF(k) — A(0)(cos k x a — cos k y a) ^^(cos k x a + cos k v a) 

(1) 
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This is not a lowest order d-wave symmetry as is usu- 
ally considered in the literature. It was Lenck and Car- 
botte, who first pointed out this factcJ. They obtained 
the superconducting gap function by using BCS theory 
with the phenomenological spin susceptibility as the pair- 
ing interaction by using fast-Fourier-transform technique, 
without any prior assumption about the symmetry of the 
gap function. They concluded that the gap structure 
although have nodal lines along k x = k y cannot have 
the simple form of cosA^a — cosk y a with a the lattice 
parameter. In a weak coupling theory language, in or- 
der to get a gap symmetry as A$F(k) one needs a pair 
potential which also has the same symmetry. We show 
in a tight binding picture by considering higher neigh- 
bours attraction that such potential is derivable up to 
third order term in equation (1). After obtaining the 
pair potential we calculate the explicit structure of gap 
function thus obtained and found to be similar to that of 
Lenck and Carbotte. We also point out that such higher 
anisotropic d-wave symmetry is key to understand larger 
2A(k) max /k B T c BCS ratio. 
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In the spirit of the tight binding description assuming 
that the overlap of orbitals in different unit cells is small, 
compared to the diagonal overlap values, the matrix ele- 
ment V{k, k') may be written as, 

3 

V(q) = Vge^ s =Vo +2^2 V n (cos q x na + cos q y na) 

8 n=1 > 

(2) 

where Rs = ±na locates nearest neighbours and further 
neighbours ; since we shall only be interested in the d- 
wave channel (i.e, in a square lattice we are not consid- 
ering 2nd, 5th etc. neighbour matrix elements as it gives 
rise to d xy , s xy channels). Thus we get, from the require- 
ment of singlet pairing symmetry i.e, A(/c) = A(— k), 

v(k, k!) = v; + £ VnfZfZ + v n gM> (3) 

n n 

where f k (g k ) = cos k x na =p cos k y na, VJ is the on-site 
term (the label r stands for repulsion, but could be at- 
tractive as well giving rise to isotropic s wave pairing) 
and the 3 rd term in equation (3) responsible for extended 
s-wave pairing will be omitted from further discussion. 

In deriving equations (2,3) we have taken into account 
attractions only along the x and y axis neighbours. How- 
ever, such attractive interaction between the 4 th neigh- 
bours also gives rise to unconventional d-wave pairing 
channel, the pair potential for the 4 th neighbour interac- 
tion that leads to singlet d-wave pairing may be obtained 
as, 

V(k, k') = 2V4(cos k x — cos k y )(l + 2 cos k x cos k y ) x 

(k — » k') + 2V4(cosfca; — cosk y )(2smk x siuk y ) 
x(k^k') (4) 

where V4 indicates strength of 4 th neighbour attraction. 

Thus, taking into consideration just the d-wave chan- 
nel, one obtains the anisotropic pair potential as, 

V(k, k') = Vxf k f v + 4V 3 f k f k ,g k g k , + 2V i f k f k ,[{\ + d k J 
x(l + d Ky ) + s Ky s Ky ] + l6V 6 f k f k > x 

~Vf k f k ,(l+g k g k ,+g 2 g 2 ,) (5) 

where we assumed 2Vi = 4V 3 = 4V4 = 16^6 =V (say) 
and considered only the appropriate contributing terms 
in V4 and Vq that leads to the gap structure in equation 
(1) [in deducing the 2 nd result of Eq (5)]. In equation (5) 
Vi, V3, V4, Vq represents strength of attraction between 
the I s ', 3 rd , 4 th , 6 th neighbours respectively and the mo- 
mentum form factors are d k;cy = 2 sin k x sin k y , s k;cy = 
2 cos k x cos k y , f k — cos k x — sin k x , g k = cos k x + cos k x , 
with = (f k ) 2 & g\ — (g k ) 2 - The actual approximation 
involved to get the second result of the above equation 
are, 



2V 4 A/*[(1 + d k J(l + d k ,J + s kxy s Ky ] ~ -f k f k , 

V e f k f k >(4gt - 2s kxy - 3)(4g 2 k , - 2s Ky - 3) ~ 

WVefkfk'ghl (6) 

These approximations are used just to retain the form of 
the gap structure (1), however, full form of the poten- 
tial, the first result of equation (5) will also be explicitly 
used. It turns out that the full potential leads to results 
very close to that obtained in the spin fluctuation the- 
ory, whereas the form given as equation (1) is just as an 
artefact of the approximations used in the refjj. 

We shall show now that the pair-potential in the sec- 
ond result of equation (5) can produce the gap symmetry 
of the spin fluctuation theory given in equation (1) up to 
3 rd order term i.e, A(k) = A(0)/ fc (l + g k + g 2 k ). (Note, in 
a weak coupling BCS theory one would tend to think of a 
potential, V(k, k') = VF k F k >, where F k = f k (l+g k + gl) 
should be essential to produce the above gap symme- 
try). Supposing that the pair-potential in the second 
result of (5) does produce the gap symmetry in (1) up 
to the third-order term, we insert the pair potential 
(the second result of equation (5) ) and the correspond- 
ing gap function into the BCS gap equation, A(fc) = 
E'kiEk')- 1 V k , k ,A(k')ttmh(f3E k ,/2). Then, comparing 
the fc-dependence from both sides of the gap equation one 
gets the required gap equation that produces the approx- 
imated spin fluctuation (ASF) gap symmetry is obtained 
as, 

A(0)=^(F/3)F fc 2 ,M tanh (^l). ( 7 ) 

k' 

This is exactly what the gap equation one would find 
using the pair potential V(k, k') — VF k F k / which cer- 
tainly produces the required gap symmetry A(k) = 
A(0)/fc(l + g k + g k ) (the pair vertex in (7) is only renor- 
malized to V/3). The symbols in equation (7) have their 
usual meanings with the superconducting quasiparticle 
energy is given by, E k — yj (e k — [i) 2 + A 2 (fc), where \i 
is the chemical potential which controls the band filling 
with the help of a number conserving equation. The tem- 
perature dependence of the chemical potential is taken 
care in the self-consistent numerical solutions of the gap 
equation. We use the band dispersion e k obtained from 
the angle resolved-ahotoemission experiment carried out 
by Norman et alj^J for the Bi— based cuprates. The 
dispersion gives experimentally measured hopping am- 
plitudes up to fifth neighbours which is plausible in the 
present calculation since further neighbours attractive 
pairing interaction is considered. 
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FIG. 1. Momentum dependence of the superconducting en- 
ergy gap function which has symmetry of that of the approx- 
imate spin fluctuation (ASF) theory (up to 3 rd order in (1)). 
The required attractive potential is derived within the ap- 
proximation (6) (a). Note the strong resemblence of this gap 
structure wvfet that obtained by Lenck and Carbotte (cf. fig- 
ure 1 of ref.cj). This higher anisotropic d-wave gap function 
yields a 2A(fe) maa; /fc_gT'c = 6 where the A(k) ma x is not at 
(k x , k y ) — (0, ±vr) but around (0, ±1.57). (b) Momentum de- 
pendence of the superconducting energy gap function when 
the attractive potential is derived exactly without the ap- 
proximation (6) leading to higher anisotropic d-wave (HAD) 
symmetry. Contrast, the deviation in the gap structure due 
the approximation (6) in (a). These gap structures (a, b) 
indeed-,consistent with that obtained by Lenck and Car- 
botte (cfJ_j). The exact HAD symmetry gap structure (b) 
although have similar structure to ASF symmetry, the gap 
becomes more sharply peaked in different k directions. Un- 
like the simplest d-wave its maximum occur around (0, ±1.41) 
which yields a 2A(fc) maa; /fcsT c = 8.15, the value srery much 
consistent with the original spin fluctuation modelul. The pa- 
eter V is such adjusted that both the gap functions have 
Tthe s^me bulk T c = 84 K at the band-filling p = 0.8. Un- 
doubtedly, these gap structures are quite different from the 
lowest order usual d-wave symmetry usually considered in the 
literature. 
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FIG. 2. (a) Variation of the gap amplitude (A(0) in eV, 
see equations (1), (7), (8)) as a function of temperature in 
Kelvin. Note, the amplitude A(0) for the usual d wave and 
the ASF (under approximation (6)) is almost same whereas 
the HAD (without the approximation (6)) is quite different. 
The gap opens up very fast below T c in the HAD case, (b) 
The BCS characteristic ratio A(k)max/k,BT c as a function of 
T/T c in different models. The solid curve represents result for 
the derived approximated spin fluctuation (ASF) symmetry 
within the approximation (6); the maximum of the gap opens 
up at a faster rate than the d-wave below T c . The dotted line 
corresponds to the HAD symmetry when no approximation 

(6) is included; this maximum gap has the fastest growth with 
lowering in temperature below T c . The dashed curve corre- 
sponds to a usual lowest order d-wave. Therefore, the special 
momentum anisotropic form in the spin fluctuation theory 
[Monthoux, Balatsky, Pines, Phys. Rev. Lett. 67, 3448] is 
one of the crucial ingredients for so high BCS characteristic 
ratio. (All the curves in the figure correspond to T c = 84 K 
at the band filling p = 0.8). 

Following the same principle as in deriving equation 

(7) , we get the gap equation for the higher anisotropic 
d wave (HAD) symmetry A(fc) = A(0)[/i(fc) + f 3 (k) + 
f4a(k) + fib{k) + fe(k)] using the full potential i.e, the 
first result of the equation (5) as, 

A(0)=X:(F/5)F fe 2 ,Mtanh(^l) (8) 

where F k = h{k) + f 3 (k) + f 4a (k) + f 4b (k) + f 6 (k) with 
fi(k) = (l/V2)/fc, h{k) = fk9k, Ua{k) = (1/V2)A(1 + 
*0> Mfc) = (1/V2) f k d kxy ), f 6 (k) = f k (gl - s kxy /2 - 
3/4). Now, we present our numerical results in figures 
1 <fc 2, for a fixed cut-off frequency f2 c = 500 K. The 
bulk T c in all the figures is fixed at T — 84 K for the 
band filling p = 0.8, which required a change in V in 
the superconducting gap equations (7), (8) and the same 
for the usual d-wave. In figures 1 (a,b) we present the 
fc-anisotropy of the gap function in the first Brillouin 
zone for k x , k y , which clearly produces d-wave like so- 
lution i.e, nodal lines along k x = k y directions as well 
as change in the sign of the gap function. However, the 
overall anisotropy is very different from the simplest low- 
est order d-w&ve (A(fc) = A(0)(cosfc x a — cosk y a)) form 
- but rather a higher anisotropic d-wave. Figure 1(a) 
presents the momentum anisotropy for the ASF model 
using gap equation (7) whereas the Figure 1(b) repre- 
sents the HAD symmetry using the gap equation (8). The 
HAD symmetry gap has sharper k-anisotropy than the 
ASF model although there is overall similarity between 
the two, namely, positions of maximum gap ar^ very 
close. Since the calculation of Lenck and CarbotteO does 
not assume any form of the superconducting gap func- 
tion and also does not include retardation effects, as in 
the original spin fluctuation modeB, we can therefore cer- 
tainly rely on comparing our results with those of Lenck 
and Capbotte. A close comparision of our results with 
those oli3 will conclusively demonstrate that the pair po- 



tential derived with distant neighbours attraction in the 
d-wave channel in the present model does produce the gap 
symmetry of the spin fluctuation theoryH. The present 
calculation thus may indicate that the phenomenological 
spin fluctuation theory includes longer range interaction 
which might be derivable from a generalised interaction 

In figure 1, the gap function shows more than one 
maximum (minimum) at the edges of the Brillouin zone, 
very similar to that obtained in referenceLa. The maxi- 
mum (minimum) is also displaced from the usual position 
((0, ±7r), (±7r, 0)) in simplest d-wavc (cf. figure 1 and 
its caption). Remarkably, this gap symmetry also pro- 
duces hiak value of 2A(fc) maa; /fcsT c = 6 same as that ob- 
tained ir£3 (cf. figure 2(b)). In Fig. 1(b) where we present 
the same as in Fig. 1(a) but without the approximation 
(6) the BCS characteristic ratio 2A(k) max / 'fcsT c = 8.15. 
These are values, close to typical of what is known for 
high-T c systemsH. Note, however—there may be some 
differences between the two (workliij and this one) only 
in details, because the band dispersion used in the two 
calculations are different which does effect pairing sym- 
metry. 

Having discussed in details the difference in the k- 
anisotropy of the (higher anisotropic d-wave) HAD sym- 
metry and that of the usual d-wave, and shown that such 
symmetries do reproduce the gap structure of the spin 
fluctuation theory, we now show explicitly in figures 2 
(a,b) that due to strong anisotropy such gaps have differ- 
ent thermal behavior in comparison to the usual d-wave 
which is principle cause for high value of 2A(k) max / ksTc- 
In Fig. 2(a) the amplitudes of higher anisotropic d-wave 
symmetries (with and without approximation (6)) and 
that of the simple d-wave are displayed as a function of 
temperature (T) — all of them have same T c — 84 K 
at a band-filling p = 0.8. In figure 2(b) we pick up the 
temperature dependencies of the maximum gap of the 
three d-wave symmetries same as shown in figure 2(a) 
and display the A(k) max /kBT c ratio as a function of the 
reduced temperature T/T c . The gap opens up below T c 
at the fastest rate for the HAD symmetry and at the 
slowest rate for the usual d-wave as the temperature is 
lowered. 

Finally, to summarize, we have derived a pair poten- 
tial from further neighbours attraction in the tight bind- 
ing scenario which produces gap symmetry of the phe- 
nomenlogical spin fluctuation theory. This study may 
particularly be justified from the fact that in models of 
spin fluctuation mediated d-wave superconductivity an 
increase in the antiferromagnetic correlation length oc- 
curs with underdoping. Such effect has also been real- 
ized very recently from angle resolved photoemmission 
(ARPES) experiment by a well known groupEj. One 
of their principle observations is that as the doping de- 
creases the maximum gap increases, but the slope of the 
gap near the nodes decreases. This feature, although 
consistent with d wave symmetry, cannot be fitted us- 
ing the simplest d-wave form of the gap but requires a 



4 



more generalised d-wave of the form, B(cos k x — cos k y ) + 
(1 — _B)(cos2fcj; — cos2fc y ) where B is a fitting parame- 
ter. (Needless to say, in the present work consideration 
of only first two terms of the first result in equation (5) 
i.e, 1 st and 3 rd neighbour interactions exactly reproduces 
this symmetry). This lead them to suggest importance 
of longer range interaction in the theory of <i-wave super- 
conductivity as one approaches the insulator. It is worth 
pointing out that such ARPES experiments in the un- 
derdoped regime measure the pseudo-gap rather than the 
truly superconducting gap. In the present calculation, a 
close look to the dispersion used in equations (7,8) will 
indicate that the Fermi surface (FS) is open in certain 
direction i.e, the FS is gapped due to pseudo-gap for- 
mation. However, if the pseudo-gap could be ascribed 
to fluctuation effects of the order parameter, then its 
value could be estimated by the mean-field BCS equa- 
tion, while the truly superconducting transition can be 
estimated only by calculations that include flutuation ef- 
fects. We thus created an example that there exist in 
nature pair potential (as we derive in real space) analo- 
gous to the spin fluctuation theory which is one of the 
leading potential theories in high temperature supercon- 
ductors, despite the fact that the principle philosophy of 
the spin fluctuation is different. We thus emphasized, 
the importance of inclusion of further neighbour attrac- 
tion in the usual d- wave, theories as is also realized in 
most recent experimenttj. The Cu — O systems being 
in a complicated circuit, the effect of columb repulsion 
may not be adequately treated with only on-site repul- 
sion and therefore, effective attractive potential may be 
achieved only after considering more distant neighbours 
terms. With this calculation we also emphasized the role 
of gap anisotropy in BCS gap ratio which may beJurthcr 
improved adopting a strong coupling approachEj. The 
higher anisotropic c?-wave symmetry as obtained in this 
work will certainly be consistent with experimental stud- 
ies in cuprates because of its similarity with ci-wave but 
will have better advantage of avoidiaft caulomb repulsion. 
We believe, this work along withcni3't3'li3 will provide new 
insight to the usual <i-wave theories of superconductivity. 
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